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I. INTRODUCTION 



One of the most important problems in the recent high energy physics is to analyze the 
quark-gluon plasma (QGP) produced in ultra-relativistic nuclear reactions JIJ. 

In a previous paper , we formulated a semi-phenomenological quantum transport the- 
ory for quark-gluon plasma fluid based on an operator-valued Langevin equation. Giving a 
mode spectrum and a damping as the input into this theory, we can easily calculate thermo- 
dynamical quantities and transport coefficients. Along this line of thought, we have already 
discussed the space-time evolution of the (l+l)-dimensional viscous quark-gluon plasma 
fluid with phase transition |§, the (3+l)-dimensional perfect fluid quark-gluon plasma || 
and the baryon-rich quark-gluon plasma |J . Furthermore we have also analyzed the existing 
experiments based on the theoretical results of the (3+l)-dimensional perfect fluid model 
with phase transition ||. 

In high energy nuclear collisions, all secondary particles, such as hadrons, leptons and 
photons, come out from a hot and dense matter which is produced at the earlier stage of 
the nuclear collision. We first naively expect that the final distribution of those particles 
directly reflects information of the hot and dense matter, as the particle source. As for 
hadrons, however, we can say that the final distribution is seriously masked by strong final 
interactions. For this reason, we prefer to observe leptons and/or photons directly coming 
from the particle source rather than hadrons themselves, in order to study the particle source 
in the earlier state at higher temperature and at higher density. On the other hand, we expect 
to observe a sort of phase transition between the hadron and quark-gluon plasma states in 
the course of cooling down of the matter from higher to low temperature. For this purpose, 
the lepton distribution would be the best to be observed, even though the observation must 
not be so easy from the practical point of view. Instead, if we can distinguish experimentally 
photons which are thermally emitted from the particle source at higher temperature and 
higher density, from photons emerging through the process tt° — > 27, the observation of 
these thermal photons would be meaningful. 
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First imagine that we have a simple case in which the average wavelength of photons 
be much longer than the size of the matter, even though this case is very far from real- 
istic cases. In this case we can expect to have a simple dipole radiation emitted from an 
oscillating dipole along the collision axis, and then to observe photons mainly distributing 
around the transverse direction. The distribution is to be compared with photons coming 
from the decay of neutral pions, which sharply distribute around the collision axis. This 
expectation must not be true, because the above assumption that the wavelength is much 
longer than the particle source is not justified. However, we could expect to observe a certain 
trend of expansion in the pr-distribution of high energy photons produced in high energy 
nuclear collisions. In this paper, we examine this kind of trend by making use of numerical 
simulations of hydrodynamical expansion of high temperature and high density quark-gluon 
plasma fluid. 

Since the thermal photon is considered to keep the information about the early stage of 
the hot matter produced by relativistic nuclear collisions, many theoretical analyses have 
already been done. Some groups P|-|TT| have analyzed the experimental data of CERN 
WA80 S+Au 200GeV/nucleon (preliminary) |T^,[rjj so as to fit their theoretical model to 
the thermal photon emission data. Most of these papers except Ref. however, dealt 



only with the photon spectrum leaving the hadron spectrum not analyzed. In this paper 
we analyze the photon and the hadron distribution produced by the hot QCD matter in a 
consistent way. 

1) We first choose parameters in the hydrodynamical model so as to reproduce the 
hadronic spectrum, i.e. the (pseudo-) rapidity distribution and the transverse momentum 
distribution. 

2) We derive the thermal production rate of photons from a unit space-time volume 
based on the finite temperature field theory. 

3) Accumulating the thermal production rate over the whole space-time region covered 
with the particle source, which is estimated by the hydrodynamical model, we evaluate the 
thermal photon distribution which is to be compared with the experimental data. 



Assuming local equilibrium for the hot QCD matter which will be produced in high 
energy nuclear collisions, we will apply the hydrodynamical models ||. As for the equa- 
tion of state, we discuss two different types of the model, i.e., the QGP fluid model with 
phase transition between QGP and hadrons and the hot hadron gas model without phase 
transition. Based on the assumption that the QGP fluid of the dominant mode obeys the 
operator- valued Langevin equation ||, we can easily deal with the thermal photon and ob- 
tain the formula for the production rate. By using the formula, we derive the thermal photon 
distribution in high-energy nuclear reactions, and compare the theoretical results given by 
the QGP fluid model with phase transition with the hot hadron gas model without phase 
transition in detail. 

In Sec. E, we shortly review the relativistic hydrodynamical model with phase transition. 
In Sec. IE, using the quantum Langevin equation, we derive the production rate of thermal 
photons from QCD matter at the high temperature region produced by high energy nuclear 
reactions. The transverse photon distribution and the asymptotic slope parameter will be 
obtained in Sec. IV. Section V is devoted to concluding remarks. 

We use the natural unit (h = c = 1 together with ks — 1) throughout this paper. 

II. HYDRODYNAMICAL MODEL WITH PHASE TRANSITION AND 

PARTICLE PRODUCTION 

The hydrodynamical equation for perfect fluid is given by 

= 0, (1) 
= (E + P)U»U U - Pg" u , (2) 

where E, P, and are, respectively, energy density, pressure, and local four velocity. 
Energy density, pressure, and entropy density are given by 

I poo 

E(T) = — p 2 dpe(p)n(e(p),T), (3) 
2ir z Jo 
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P(T) = ± r p 2 dp^-Me(p),T), (4) 
6n 2 Jo e(p) 

E + P 

S(T) = -±r~. (5) 

Following Ref. M, we assume a simple model for the mode spectrum of the fluid 

1 - tanh , 1 + tanh ^=S? 



e(p) = Ayjp* + d —+ | p | (6) 

Here we suppose that the fluid in the QGP phase is dominantly composed by u-, d-, s-quarks 
and gluons and that the fluid in the hadron phase is dominantly composed by pions and 
kaons. In this case we put A = 1.89, M = 200 MeV, T c = 160 MeV, and d = 2 MeV based 
on a previous analysis. With these parameters, we obtain the phase transition-like behavior 
of energy density (see Fig. 1), which seems to reproduce the Lattice QCD result 0. 

In a previous paper ||, by making use of the following two models: 1) the QGP fluid 
model with phase transition, 2) the hot hadron gas model without phase transition, we have 
analyzed the pseudo-rapidity distribution of charged hadrons in S+Au 200 GeV/nucleon 



collision obtained by CERN WA80 |JL2|| . In this paper in order to develop the hydrodynam- 
ical model furthermore, we are going to analyze the p-p-distribution of neutral pions also 
given by CERN WA80 |T3| as well as the pseudo-rapidity distribution of charged hadrons. 



According to the previous analysis ||, we use the first model (the QGP fluid model with 
phase transition) specified by the initial temperature Tj = 195 MeV, the critical temperature 
T c = 160 MeV, and the freeze-out temperature Tf = 140 MeV, and the second model (the 
hot hadron gas model without phase transition) specified by Tj = 400 MeV and Tf = 140 
MeV. We choose other parameters being the same as given by Ref. 0. For these models, 
we obtain theoretical results of the hadronic spectrum. See Fig. 2 and Fig. 3. From Fig. 2 
and Fig. 3, we observe that the both models can consistently reproduce the experimental 
data. The hot hadron gas model should have the initial energy of the fluid exceeding the 
total collision energy 0, so that we conclude that the hadron gas model is not accepted. 
We discuss the thermal photon distribution based on these two models in Section IV, in 
which we will see that the hot hadron gas model will fail again in the analyses of the photon 



experimental data. 



III. THERMAL PRODUCTION RATE OF PHOTONS 

Transition amplitude for a process involving emission of a photon with momentum k and 
polarization from the QCD matter is given by 

T Pa = (/3;fc,£ (A) ;out | a;0;in), (7) 

where | a; 0; in) and | /3; k, e^; out) stand for, respectively, the initial state and the final 
state. Usual technique of the reduction formula enables us to give the following transition 
amplitude 

Tpa = (A out | a { X J t (k) | a; in) 

= ig>*eW I e ikx ((3; out | f H (x) \ a; in) 

J ^{2^2k 

= ig XX 'e^ I . e*-(fl in | T{Sj?{z)) \ a; in), (8) 

J yj(27t) 3 2k 

from which we obtain transition probability per unit space-time volume 
R{k,e^) = ^ f \T Pa | 2 



y (27r) 3 2fc 

x (a; in | f(tftf(x)) | P;in)((3;in | T(Sj!(y)) \ a- in). (9) 

Since we are interested in one photon state, we should sum up with respect to the final 
states of QCD matter denoted by (3. To the lowest order in QED, we obtain 

RfaeM) = A^- / ^^f lk{X - V) ^™ I r(x)f(y) I a; in). (10) 

Assuming that a certain mode is dominantly excited in a local equilibrium system of 
the hot QCD matter and the canonical operator of that mode obeys the quantum Langevin 
equation (see Appendix), then we can replace the above matrix element with the ensemble 
average in the sense of quantum Langevin equation, i.e., 
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(a; in | f{x)j u {y) \ a; in) =>• (jQ. L .{x)j^ L {y)) Q . L ., (11) 

where subscript Q.L. represents the ensemble average. Through this paper, we will omit the 
subscript Q.L. for the simplicity. 

In the case of bosonic mode (pion or kaon) in the hadron phase, the source function is 
given by 

f(x)=i:<P^x)(B 1 -d> x )</)(x):. (12) 
The boson field operator in the theory of the quantum Langevin equation is represented as 

4>{ X ) = [ . d3p r dw\P^- P {p,u) 

< A(p,u) B\p,u) ■ \ 

X [uj-E( P ,T) e + u J -E*(p,Tf J' (13) 

E(p,T)=e(p)-^r(T), (14) 

where e(p) = \/p l + M' 2 and j(T) = cT are, respectively, the mode spectrum and the 
damping, which we have given to the theory as input. The free parameter c in j(T) is so 
chosen as to fit the photon transverse momentum distribution to the experimental data: we 
have determined c = 0.01. In the case of fermionic mode (quark) in the QGP phase, a source 
function is given by 

f{x) =: $(x)f-il)(x) :, (15) 

where 



Here we have put e(p) — \ p \ and 7(T) = cT, with c = 0.01. 

Summing up with respect to polarization and wave number, we can obtain the production 
rate 



R(T)= [ d 3 kJ2R(k,£ W ), (17) 

J A 

which means the number of the emitted photons per unit space-time volume of the hot 
QCD matter. In order to evaluate the production rate at temperature T, let us introduce 
the following simple model 

1 - tanh 1 + tanh 
R(T) = Rhad h Rqgp ^ , (18) 

where we have supposed that the phase transition takes place around critical temperature 
T c within width d. In the infinite T c limit, we can only keep the hot hadron gas model. 
Based on the Langevin technique mentioned in Appendix, we can easily calculate Rhad and 
Rqgp as 

d 3 k r d 3 p 1 r°° dwi 7(T) 



Rhad " E e ' / (2vr )3 2 A; / (2% I 



(2tt) 3 2A;o J (2tt) 4 Jo 2uj x \ u x - E(p 1 , T) 



2 



X 



d 3 p 2 f°°duj 2 7 (T) 



(27r) 4 y 2lu 2 I oo 2 - E(p 2: T) 



2 



x (27r) 4 (-n 6 (^){l + n b (uj 2 )}( Pl + p 2 ) 2 5 4 ( Pl - p 2 - k) 
-{1 + n b (uj 1 )}n b (u 2 )(p 1 + p 2 ) 2 5 4 (p 1 -p 2 + k) 
- n b (ui) n b (uj 2 ) (pi -p 2 ) 2 5 4 (pi +p 2 - k)), (19) 



d v o 2 f d 3 k r d 3 Pl r°° du x 7 (T) 

n QGP ^ 6e s J ( 2 7t) 3 2A; J (2vr) 4 Jo 2u ± \ u x - E( Pl , T) 

d 3 p 2 f°° du 2 7(T) 



2 



(2tt) 4 7 2uj 2 I uj 2 - E(p 2 ,T) | 2 
x(2 7 r) 4 (n / (cc; 1 ){l - n f (io 2 )}(8 Pl p 2 - 16^^)5^ -p 2 -k) 

+{1 - n/(wi)}n/(w 2 )(8pip 2 - ^6\fl%\fpll)S 4 {Pi ~ P2 + k) 

+ rifM n f (u 2 ) (8 Pl p 2 + 16 Jp 1 1 Ji%)5 4 (p 1 +p 2 -k)). (20) 



Here the summation s is over the source particle, and e s stands for the electric charge of a 
source particle. These formulas contain three kinds of photon emission processes, as shown 
in Fig. 4. The double line stands for off-shell particles in thermal bath. Taking into account 
the structure of Rhad and Rqgp having factors r y(T)/ \ uj — E(p,T) | 2 , we can understand 



that the particle spectrum in heat bath is spread due to the presence of random force. The 
production rate as a function of temperature T is shown in Fig. 5. 

In order to obtain the theoretical formula for hydrodynamical quantities and transport- 
theoretical coefficient, we have assumed that the local system is in uniform thermal- 
equilibrium on a microscopic space-time scale. On the other hand, the QGP fluid hydrody- 
namically evolves in the macroscopic space-time region, specified by space-time coordinate 
x. Solving the hydrodynamical equation, we obtain the x-dependence temperature T(x) 
and local four velocity U IJ '(x). Following the well-known procedure, we first relate the pro- 
duction rate in each local system of the QGP fluid with one in the whole center-of-mass 
system of the QGP fluid. Since the number of photons emitted from the local system is 
Lorentz- invariant, we obtain the following relation between the local system at space-time 
point x and the overall center-of-mass system 



d 3 R c . m . _ d 3 R(T(x)) 
dk 3 dk' 3 



(21) 

k' Q =U^{x)k„ 



Integrating the above R c . m . over the whole space-time volume in which the particle source 
exists, we obtain momentum and transverse momentum distributions 



d 3 N _ f d ^k' d 3 R(T(x)) 



dk 3 J k dk 



/3 



(22) 

k' =U^(x)k„ 



1 dN r „ d 3 N /7 , JN 

V T W T = J d^ V (k T ,k L ,cf>), (23) 

which are to be compared with experimental data. Here temperature T(x) and local four ve- 
locity U^{x) at space-time point x are given by the numerical solution of the hydrodynamical 
model. 



IV. RESULTS AND DISCUSSION 



Figure 6 shows the numerical results of Eq. (23) compared with the experimental data 
(S+Au 200 GeV/nucleon collision) obtained by CERN WA80 flyj. The solid curve and the 
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dotted curve are, respectively, the contribution of the hadron phase region and that of the 
QGP phase region in the QGP fluid model. The whole thermal photon distribution given by 
our QGP fluid (with phase transition) model is the sum of them, but we can easily see that 
the contribution of the QGP phase region is negligibly small. Using our hydrodynamical 
model for the QGP fluid with phase transition, we also see in Fig. 6 that we reproduce the 
experimental data of WA80 well. The dashed curve stands for the photon distribution given 
by our hot hadron gas model. The theoretical curve (the dashed curve) deviates from the 
experimental data in both absolute value and slope. Our formalism keeps a free parameter 
c in the damping function which is directly related to the intensity of the photon emission 
but which is almost independent of the slope. Even if we choose another value for a free 
parameter c, the hot hadron gas model cannot reproduce the experimental data due to the 
deviation of the slope. Therefore, we say that we reproduce the WA80 experimental data 
consistently only with the QGP fluid model with phase transition. 

Needless to say, the initial temperature or the phase transition temperature are very 
helpful to understand high energy nuclear reactions, if they are derived from theoretical 
analyses as mentioned above. We naively expect the thermal photon distribution to directly 
reflect the characteristic temperature, but we have to pay attention to the following points: 
1) temperature T cannot be invariant under Lorentz transformations, as was mentioned in 
Sec. II, 2) each local system in the QGP fluid has different temperature. For these reasons, we 
introduce a little mathematical manipulation in the following way. The thermal distribution 
from the volume element with velocity v and temperature T should read as 



expi-- T 



=>• ex p( — 

local system 



cm. system 



k U - k L U L - k T U R cos9 
= exp( ), (24) 

where = (-y=y, -y=y)- In order to pick up the most dominant contribution to the 
transverse momentum distribution, we put fcx = and rewrite Eq. (24) as 

exp( ^ ) = exp( j = — —), (25) 

1 V L / 1+"T T 
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by which we can define the effective temperature T e ff of the fluid at the volume element 
with velocity vl,vt and temperature T by 



Furthermore we have assumed that the largest value of T e // dominates in the slope of the 
transverse momentum distribution for extremely large kx region. The transverse momentum 
distribution at rapidity y = and its slope parameters are, respectively, shown in Fig. 7 and 
Fig. 8. Table 1 shows the maximum T e ff given by our numerical results of hydrodynamical 
simulation and the slope parameter at kx = 20 GeV for the above two models. We see in 
Fig. 8 that the slope parameter T s of each model asymptotically tends to maximum value 
of T e ff shown in Table 1. Through comparison of T e ff evaluated by the numerical results 
of hydrodynamical simulation with the asymptotic slope parameter T s of the transverse 
momentum distribution in Table 1, we know that the asymptotic slope parameter T s has 
possible origins different from each other for the above two models: The critical temperature 
dominates the asymptotic slope parameter in the QGP fluid model with phase transition, 
while the initial temperature dominates the asymptotic slope parameter in the hot hadron 
gas model without phase transition. 



We have derived the thermal photon distribution emitted from a hot matter produced 
by the high energy nuclear collisions, based on hydrodynamical model, and compared these 
theoretical results with S+Au 200 GeV/nucleon data obtained by CERN WA80. We have 
observed that only the QGP fluid (with phase transition) model can consistently reproduce 
the above experimental data. Furthermore we have determined the asymptotic slope pa- 
rameter of the transverse photon distribution through comparison of our numerical results 
of hydrodynamical simulation with the thermal photon distribution. 




(26) 



V. CONCLUDING REMARKS 
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APPENDIX : QUANTUM LANGEVIN EQUATION 

We introduce a quantum Langevin equation for the annihilation operator 

ij t a{p, t) = f K{p, t>(p, t')dt' + f{p, t), (27) 

and its hermit conjugate, where K(p,t) and f(p,t) are input parameter and random force, 
respectively. Canonical operators a(p, t) and a)(p,t) should satisfy following two require- 
ments: 

1) Equal time commutation relation 

[a(p,t),J{p,t)] T = 6 3 (p-p'), (28) 

where =F stand for boson and fermion respectively. 

2) The Kubo-Martin-Schwinger condition 

(a\p,t)a(p',t' + z(3)) = (a(p',t')a\p,t))e^. (29) 

We can fix the detailed properties of a and / so as to satisfy above requirements. 
The random force operator can be expanded as 

f(p,t) = /^d^[^^] 1 ^(p ) aOA(p ) aOe-* y *, (30) 

where A(p, uo) is a canonical random operator whose canonical commutation relation is given 
by 

[A(p,u),A^p',J)] =6 3 (p-p')6(u-u;'). (31) 
The stationary solution of Eq. ([27]) is given as 
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a( P , t) = ^ MP^P(P.-)MP^ e -^ (32) 
J-oo V 27T uJ — E(p,u) 

The Kubo-Martin-Schwinger condition (p9|) can be satisfied by the following ansatz: 



(A(p,u)) = (i4t(p,w)) = l 

(A(p,c)At(p',c')) = 5 3 (p-po%-^)[i+e^,r)], 

(A{p,u)A(p,u)) = (A^(p,uj)A^p,uj)) = 0, 



(33) 



where £ takes +1 for boson and —1 for fermion. 

From (|30| ) and ([33]), quantum fluctuation-dissipation theorem for the operator- valued 
Langevin equation (p7[) is obtained 



(/(fc,f)> = (/t(p,t)> = 

(P(P,t)f(p',t>)) = 5 3 (p-p')r^ 2 (P^)T(P^)n(a;,r)e-^-*') 

(/(P. *)/V. *')> = ^ 3 (P " P') jf^fo w fr(P, w) [1 + £n(w, T)]e +i ^ 

(f(p,t)f( P ',t')) = (p(p,t)fKp',t')) = o 



We have formulated a quantum Langevin equation consistently. We make thermal field 
operators for bosons (|T3| ) and fermions ([16]) by means of the above operator a(p, t) and 
a! fat). 
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Table 1 



MdHpI 


T (MeV") 

J. 1 1VJ.C V 1 


Urn 


1) T 


T tt (Me VI 


T fMpV") 

± s ^IVIC V ) 


QGP fluid (QGP phase) 


157.5 


0.53 


0.11 


280.8 


273.2 


QGP fluid (hadron phase) 


157.5 


0.53 


0.11 


280.8 


270.5 


Hot hadron gas 


400.0 


0.0 


0.0 


400.0 


390.1 



TABLE CAPTION 

The maximum T e ff for the QGP fluid model and the hot hadron gas model in our 
hydrodynamical simulation, and the slope parameters T s at kx = 20 GeV in Fig. 8. T e ff is 
evaluated by Eq. (26). 

FIGURE CAPTION 

FIG. 1 The phase transition-like behavior of energy density distribution. The solid curve 
stands for our QGP fluid model with phase transition, the dotted curve for the fluid of the 
hot hadron gas composed of massive pions and kaons, and the dashed curve for the fluid 
of the quark-gluon plasma composed of massless u-, d- s-quarks and gluons. The critical 
temperature T c = 160 MeV. 

FIG. 2 The pseudo-rapidity distribution of charged hadrons in S+Au 200 GeV/nucleon 
collision. The experimental data was obtained by CERN WA80[12]. The solid curve and 
the dashed curve stand for, respectively, the QGP fluid model with phase transition and the 
hot hadron gas model without phase transition. 

FIG. 3 The transverse momentum distribution of neutral pions in S+Au 200 GeV/nucleon 
collision. The experimental data was obtained by CERN WA80[13]. The solid curve and 
the dashed curve stand for, respectively, the phase transition model and the hot hadron gas 
model. 

FIG. 4 Three kinds of photon emission processes represented by Eqs. (19) and (20). The 
double line stands for off-shell particles in the heat bath. The wavy line stands for the 
external line of the photon. 
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FIG. 5 The production rate as a function of temperature represented by Eq. (18). The 
solid curve stands for the phase transition model, the dotted curve for the production rate 
of the hadron phase and the dashed curve for the production rate of quark-gluon plasma 
phase. The critical temperature T c = 160 MeV. 

FIG. 6 The transverse momentum distribution of thermal photon in S+Au 200 
GeV/nucleon collision. The upper limit of experimental data was estimated by CERN 
WA80[15]. The solid curve and the dotted curve are, respectively, the contribution of the 
hadron phase region and the QGP phase region in the QGP fluid model with phase tran- 
sition. The dashed curve is the photon distribution of the hot hadron gas model without 
phase transition. 

FIG. 7 The transverse momentum distribution at rapidity y — 0. The solid curve and 
the dotted curve are, respectively, the contribution of the hadron phase region and the QGP 
phase region in the phase transition model. The dashed curve is the photon distribution of 
the hot hadron gas model. 

FIG. 8 The slope parameter of the transverse momentum distribution which is shown 
in Fig. 7. The solid curve and the dotted curve stand for, respectively, the hadron phase 
and the QGP phase in the QGP fluid model with phase transition. The dashed curve is 
the photon distribution of the hot hadron gas model without phase transition. The solid 
bold line and the dashed bold line stand for, respectively, T e ff of the QGP fluid model and 
T e ff of the hot hadron gas model. The values of the asymptotic slope parameter T e ff are 
evaluated by Eq. (26) and shown in Table 1. 
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Fig. 2 
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